ds = read.csv("http://www.math.smith.edu/sasr/datasets/bls.csv",header=FALSE)
jobs = with(ds, cbind(V2, V3, V4, V5, V6, V7, V8, V9, V10,V11, V12, V13))
jobsts = ts(as.vector(t(jobs)), start=c(1945, 1), frequency=12)
plot(jobsts, plot.type="single", col=4,ylab="number of jobs (in thousands)")

abline(h=0)
presline = function(date,line,name)
{
  mtext(at=date,text= name, line=line)
  abline(v = date)
}

presline(1946,1,"Truman")
presline(1953,2,"Eisenhower")
presline(1961,1,"Kennedy/Johnson")
presline(1969,2,"Nixon/Ford")
presline(1977,1,"Carter")
presline(1981,2,"Reagan")
presline(1989,1,"GHW Bush")
presline(1993,2,"Clinton")
presline(2001,1,"GW Bush")
presline(2009,2,"Obama")
